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Monte Carlo simulation of the charge carrier transport in disordered nonpolar organic materials has been 
canied out. As a suitable model we considered the model of quadrupolar glass. A general formula for the 
temperature and field dependence of the mobility was suggested. A comparison with experimental data has 
been carried out. 
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1 Introduction 

Characteristic feature of disordered organic materials is a long range spatial correlation in the distribution 
of random energies U{r) of transport sites IT] ID. This particular feature stems from two fundamental 
properties of typical organic glasses: first, they have very low concentration of intrinsic free carriers, and, 
second, organic glasses usually have high concentration of dipoles and quadrupoles. Due to almost zero 
concentration of free carriers and lack of screening such molecules provide long range contribution to 
the overall energetic disorder for charge carriers. Long range sources inevitably lead to the strong spatial 
correlation of the resulting random energy landscape. Experimental data and realistic models of organic 
glasses suggest that the density of states usually has a Gaussian form with typical rms disorder a ~ 0.1 
eV ISjIHISl. Hence, the only relevant characteristic of the random energy U{f) is the binary correlation 
function C{f) = (J7(r)J7(0)). Exact solution of the ID transport model shows that this function directly 
dictates a functional form of the drift mobility field dependence ij.{E) For this reason organic glasses 
with different C{r) must have different mobility field dependences. 
For example, if correlation function decays as a power law 



C(f) = Aa' j (1) 

(here a is a minimal distance between transport sites), then the ID transport model gives in the most 
important case of strong disorder (7/3 ^ 1 Q 

/ 1\ 1 /eai?\"^ 

\nfi/ Ho ^ ~'2aa - (3^ + I l + a /3{Ana 13)^ i—j , ^ = l/fcT, (2) 

where a is a radius of the decay of the wave function of the transport level. This particular case is of the ut- 
most importance, because in the simplest model of polar disordered organic materials (the model of dipolar 
glass lUl) the correlation function decays as l/r, and in the simplest model of nonpolar disordered organic 
materials (the model of quadrupolar glass (QG) Q) C(r) oc Simplest realizations of the models 
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could be considered as regular lattices with sites occupied by randomly oriented dipoles or quadrupoles, 
respectively. 

According to the result of the ID model, in polar organic materials the mobihty field dependence should 
obey the law In /i oc iJ^/^, and in nonpolar materials the dependence is In oc E"^/^. This result shows 
that the earlier belief, promoted by the Gaussian Disorder Model (GDM) [Sl, that there is a universal 
mobility field dependence in disordered organic materials is, in fact, wrong. The GDM assumes that there 
is no spatial correlation in the random energy landscape and corresponding mobility field and temperature 
dependence could be formally considered as the particular limiting case of (|2]i for rt — > cxd. 3D Monte 
Carlo simulation generally supports the result of the simplified ID model and typically only modifies the 
numerical coefficients in Q fSl. 

At the same time, experiment indicates that the universal Poole-Frenkel (PF) mobility field dependence 
In /i oc E^/"^ is usually a good approximation for the description of the time-of-flight (TOP) data. The 
major goal of this paper is to present results of the extensive 3D Monte Carlo simulation for the QG 
model and discuss how well they fit the TOP data for charge transport in nonpolar organic materials, and 
how the intrinsic mobility field dependence In /i oc i?'^/^ can be reconciled with the seemingly universal 
Poole-Prenkel dependence In /i oc i?^/^, found in experiments. 



2 Quadrupolar glass model: transport simulations and results 

Simulations have been done for a lattice QG model with the size of basic sample of 128 x 128 x 128 
sites of a simple cubic lattice with the lattice scale a and periodic boundary conditions imposed. Particular 
distributions of U (r) have been generated in the usual way. There is no correlation among the fluctuations 
in momentum space U{k) for different fc, so we generated random Gaussian field U{k) for the quadrupolar 
correlated function, which with very good accuracy can be described by ([U with n = 3 and A w 0.5 for 
r > a, and, of course, C(0) = cr^ Q, and then calculated Pourier transform to get U{r). Carrier hopping 
has been simulated using the Miller- Abrahams hopping rate 




which is believed to be a good approximation for the hopping process in organic glasses ||3]; we used 
aa = 5, as in |[3][11. In OJ random energies Ui include the shift from the applied field E. 

Details of the Monte Carlo simulation and analysis of the transport data are almost identical to those 
described in fSl. General features of the current transients are usual: if cr/? is not very large and the transport 
layer not too thin, then a well defined plateau is developed, demonstrating essentially quasi-equilibrium 
charge transport with average carrier velocity independent of time. In the quasi-equilibrium regime the 
dependence of /i on T and E is shown in Pig. [T] It was found that the general field dependence of the 
quasi-equilibrium mobility closely follows prediction of ID model (i.e.. In /x oc E^^^), though the numeric 
coefficients are different. A phenomenological relation of /i on T and E could be described as 

/i = /^o exp (-2aa) exp -O.STct^ + Cq {d^'^ - Pq) {eaE/af^ , a = a/kT, (4) 

with Cq « 0.87 and Pg ?a 1.91. 



3 Quadrupolar glass: comparison with experimental data 

Close inspection of Pig. explains, why the description of the mobility field dependence in nonpolar 
materials in the usual PP presentation In vs. E'^/^ provides a good hnear fit in many cases. Indeed, if 
field range is not too wide (one order of magnitude or even more narrow, which is typical for such materials 
[9 10 11 J2J ), then it is nearly impossible to distinguish the true quadrupolar dependence (3) from the PP 
one (see Pig. |2]l. 
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Fig. 1 Mobility field dependence in the QG model for different values of kT /a (from the top curve downward); 
straight lines in (a) indicate the fit for (|4); ^J^s = /io exp (— 2aa). Plot of the simulation data in the usual PF presentation 
In /i vs E^^'^ (b) demonstrates much stronger deviation from the linearity in the weak field region (straight lines serve 
as a guide for an eye). If a « 1 nm and a ~ 0.1 eV, then eaE/a ~ 1 for _E = 1 x 10^ V/cm. 



In this situation the use of the proper method of the mobility calculation from the TOF curve is extremely 
important. There are two most frequently used procedures for the calculation of ^ = L/Et from the TOF 
data in double linear current versus time plot: 1) calculation that uses time ti determined by the intersection 
of asymptotes to the plateau and trailing edge of the transient and 2) calculation that uses time ti/2 for 
current to reach half of its plateau value; here L is a thickness of the transport layer Unfortunately, the first 
method is a method of choice for most papers. One can find in literature the statement that the difference 
between two procedures is not very important for determination of the temperature and field dependence 
of the mobility. This is not true. Sometimes the first procedure distorts the functional type of the mobility 
field dependence lfT3l[T4l . Procedure that uses ti/2 is much more reliable and generally agrees well with 
the standard definition of the mobility as /i = v/E, where v is an average carrier velocity llT4l . This very 
procedure must be used for the reliable estimation of the mobility field dependence. 

Let us discuss how good is the model of quadrupolar glass for the description of the transport properties 
of nonpolar organic materials. First of all, most transport molecules, polymer binders or inert dopants do 
possess dipole moments, albeit sometimes quite small ones. For example, in a typical nonpolar polymer 
matrix, polysterene, a unit of the polymer chain has a small dipole moment 0.4 D IfTSl . Hence, a more 
proper model for many nonpolar glasses should be a mixture of the dipolar and quadrupolar glasses. Un- 
fortunately, such model should be hardly tractable because the mobility field dependence is a mixture of 
dipolar and quadrupolar contributions and individual contributions cannot be easily separated even in the 
ID case lfT6l . 

Second, there is another, more fundamental limitation for the direct applicability of the QG model. In 
glasses there is no translational order and, as a consequence, the symmetry is broken. For this reason 
molecules could possess induced dipole moments; for example, in symmetric PPVs, having no permanent 
dipole moment, induced dipole moments could reach ~ 0.5 D ifTTl . Induced dipolar disorder is usually 
smaller than the quadrupolar disorder, and yet for the weak field region its contribution can be quite com- 
parable to the contribution of the dominant quadrupolar disorder ||6l. In this situation no simple relation 
for the mobility field dependence can be derived. Very probably, the only tool suitable to provide the mo- 
bility field and temperature dependence in this complicated situation is a direct simulation of the materials 
structure with subsequent calculation of all site energies and transfer integrals ifTSi . 
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Fig. 2 Mobility field dependence in polysterene doped with 50% of weakly polar transport dopant TASB (bis(ditolyl- 
aminostyryl)benzene, dipole moment 0.54 D), plotted as In /i vs E^^"^ in accordance with (|4]l and as a usual In ^ vs 
^1/2 pp pjpj Straight lines show the best linear fits and numbers indicate the temperature. Both plots demonstrate 
approximately the same linearity. The experimental data has been taken from Figure 3 in II12I . 

4 Conclusion 

We have shown that the field dependence of the mobility in nonpolar organic materials differs from the 
corresponding dependence in polar organic materials (Poole-Frenkel dependence). Yet the difference is not 
very large and this is the reason why the PF dependence usually fits well experimental data for nonpolar 
materials too. Nonetheless, evaluation of the disorder parameters (i.e., a) for nonpolar materials from the 
PF mobility dependence can give inaccurate results. 
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